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ABSTRACT 

We consider the possibihty that sohd molecular hydrogen is present in interstellar space. If so 
cosmic-rays and energetic photons cause ionisation in the solid leading to the formation of . This 
ion is not produced by gas-phase reactions and its radiative transitions therefore provide a signature of 
solid H2 in the astrophysical context. The vibrational transitions of are yet to be observed in the 
laboratory, but we have characterised them in a quantum-theoretical treatment of the molecule; our 
calculations include anharmonic corrections, which are large. Here we report on those calculations 
and compare our results with astronomical data. In addition to the H^ isotopomer, we focus on 
the deuterated species (HD)j]' which is expected to dominate at low ionisation rates as a result of 
isotopic condensation reactions. We can reliably predict the frequencies of the fundamental bands 
for five modes of vibration. For (HD)j^ all of these are found to lie close to some of the strongest 
of the pervasive mid-infrared astronomical emission bands, making it difficult to exclude hydrogen 
precipitates on observational grounds. By the same token these results suggest that (HD)j|" could 
be the carrier of the observed bands. We consider this possibility within the broader picture of 
ISM photo-processes and we conclude that solid hydrogen may indeed be abundant in astrophysical 
environments. 

Subject headings: ISM: lines and bands — line: identification — molecular processes — dust, extinction 



1. INTRODUCTION 



Astronomers observe matter in a great variety of phys- 
ical conditions and, correspondingly, hydrogen exists in a 
variety of states in different astronomical contexts. Liq- 
uid metallic hydrogen is inferred to be present i n a degen- 



erate core in th e planets Jupiter and Saturn (Hubbard, 



Burrows and Lunine,2002j , and in the interstellar context 
there are known to be three phases of hydrogen in rough 
pressure equilibrium: hot, ionised, low-density plasma; 
atomic gas; and molecular hydrogen in cold, dense, in- 
terstellar clouds. Might there also be solid molecular 
hydrogen? 

That question was first addressed in the late 1960's 
when it was suggested that thick mantles of solid H2 
might be deposite d onto the surfaces of graphite grains — 
see, for exa mple, Wickramasinghe and Reddish fl968j); 
Hoyle , Wickram asingh e and Reddish (1968); Wickranias^ 
mghe and Krishna Swamy (1969f~T?ifrTShese early pa- 
pers employed a poor description of the H2 saturated 
vapour pressure and in fact solid H2 mantle formation is 
not expected under the conditions known to exist in th e 



diffuse ISM (Greenberg and de Jong 1969 Field 1969). 



The difficulty is simply that the pressure in the diffuse 
interstellar medium is much smaller than the saturated 
vapour pressure of H2, even for temperatures as low as 
that of the microwave background. 

The argument just given may, however, be misleading 
because the phase equilibrium of H2 in astrophysical en- 
vironments may differ significantly from that of pure H2 . 



In particular the electrostatic binding energy contributed 
by ions in the lattice may substantially increase the total 
binding energy of the solid, even at very low ion concen- 
trations. Astrophysical H2-ice may thus be more robust 
than is commonly supposed. 

It also remains a possibility that there exists gas which 
is so cold and dense that pure hydrogen precipitates can 
form. But that gas cannot be part of the diffuse ISM 
and must, instead, form distinct, self-gravitating entities. 
Clouds of this type - cold, dense, mostly transparent and 
having very little thermal emission - were postulated by 
Pfenniger and Combes (1994), Pfenniger, Combes and 
Martinet (1994), to be a significant component of the 
Galaxy's dark matter. That picture was subsequently 
considered by many other authors (e.g. Gerhard and 
Silk 1996, Walker and Wardle 1998, Draine 1998, Sciama 
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2000, Walker 2007). It was argued by |Wardle and Walker 
(|1999|) that small amounts of solid H2 can confer thermal 
stability on these clouds. 

Irrespective of these theoretical deliberations, obser- 
vation provides the final arbiter of whether solid H2 is 
present in the Galaxy. To that end we must first ask 
how solid H2 might be detected if it were present? That 
is the question addressed in this paper. It is not a triv- 
ial matter because solid H2, like the gas from which it 
forms, is cold, dense and very nearly transparent in the 
optical. We make no attempt at a comprehensive assess- 
ment of the detection problem; instead we focus on the 
properties of a particular molecule, H^, which is known 
to form in solid H2 exposed to ionising radiation (see 
§2). Cosmic-rays and energetic photons pervade inter- 
stellar space and we therefore expect H^ to be present 
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at some level in solid, interstellar H2. Furthermore, this 
molecular ion does not form in any significant quantities 
in the gas-phase chemistry of H2, so it is a good tracer 
of solid H2 . 

Currently H^j" is not well characterised experimentally 
so in order to predict its astrophysical manifestations we 
are first obliged to confront problems in quantum chem- 
istry. To tackle these we have utilised the Q-Chenj^ soft- 
ware package (Shao et al. 2006), with some bespoke 
extensions, for our calculations. For the most part we 
have concentrated on the bare H^i" molecule, whereas in 
solid H2 it is expected that the ion forms many weak 
bonds with surrounding H2 molecules, leading to a much 
larger complex. Our focus on bare is motivated by 
the fact that solvati on does not greatlv disturb the stru c- 
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ture of the H, 

and thus the high frequency vibrational transitions of the 
solvated forms should closely resemble those of the bare 
molecule. Moreover we do not yet know exactly what 
form the solvation takes. We do, however, give some 
consideration to how the properties of the solvated forms 
differ from those of the bare molecule, and where possi- 
ble we estimate the size of these effects. This has been 
done by reference to the properties of 11+4 — Hfj" (112)4 
— a cluster which is much smaller than those which are 
likely to occur in practice. 

We have neglected any other interactions between H^i" 
and the host lattice. This is a reasonable assumption 
for an initial exploration as the intermolecular forces in 
solid hydrogen are very feeble. Indeed solid para-H2 is 
commonly used for matrix isolation spectroscopy of small 
molecules (see, e.g.. Oka 1993, Anderson et al 2002), and 
this is precisely because it permits nearly free vibration 
and rotation of the impurity species. In the astronomical 
context there is plenty of time for ortho-para conversion 
to take place, so any cold H2 is expected to be predom- 
inantly in the rotational ground state and is thus in the 
"para" form. 

This paper is structured as follows. In the next section 
we discuss relevant laboratory results on irradiated solid 
hydrogens and their interpretation. In §3 we present our 
theoretical model of the Hr 
^+ 



j-g molecule and its isotopic 
variant (IID)j]". The latter is expected to be the most 



abundant isotopomer at low ionisation rates whereas the 
former dominates at high ionisation rates. Section 4 
compares the computed vibrational transitions of 
and (HD);^ with mid-IR astronomical observations. The 
properties of ten other isotopomers, which may be of in- 
terest in laboratory studies, are tabulated in Appendix 
A. Similarities between computed and observed spectra 
prompt us to suggest that (HD);^ may actually be re- 
sponsible for the observed mid-IR astronomical emission 
bands; the implications of this suggestion are discussed 
in §5 and we offer our conclusions in §6. 

2. LAB STUDIES OF H+ IN IRRADIATED, SOLID H2 

At temperatures less than the triple-point (13.8 K), 
hydrogen is able to form a low-density solid {p ~ 
0.086 gcm""^). It is a "quantum solid", meaning that 
the zero-point motion of the molecules is a large fraction 
of the intermolecular separation. Solid hydrogens have 
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been extensively studied in the laboratory; see |Silvera| 
^980) and Souers (1986| for reviews of their properties; 
p^n Kranendonk, ((1983| provides theoretical models for 
some of these properties. 

Although the natural abundance of tritium is negligi- 
ble, its importance in fusion research led to studies of 
solid hydrogens in which tritium is incorporated. These 
solids naturally exhibit the effects of irradiation because 
tritium undergoes beta decay to ^He with a half-life of 
approximately 12 years. The energetic electrons which 
result from beta decay cause ionisation and dissociation 
of the H2. Various radiation-induced phenomena have 
received attention in the literature, as described below. 

2.1. Radiation-induced absorption lines 

In the first studies of irradiated solid hydrogens it 
was discovered that some extra lines were present in 
the absorption spectru m near the funda.menta l vibra- 
tional band of the solid (Souers et al||1980 1981). These 



lines were attributed to the presence of ions m the ma- 
trix which induce a Stark-shift in the quantum states of 
neighbouring molecules. 
Detailed studies of these spectral features were made 



by Brooks et al (11983^ 1985^ using proton-beam irradia- 
tion of tritium-tree solid U2 — an approach which allows 
complete control over the magnitude and timing of the 
radiation doses. The radiation-induced lines were seen to 
persist for many hours after irradiation ceased, i mplying 
that the ions involved have very low mobilities. |Brooks| 
et al ( 1985) identified the immobile anion as an "electron 
bubble" (i.e. an electron trapped in a lattice defect) but 
did not claim to have identified the immobile cation; they 
suggested Dj^, or the solvated form Dg = 0^(02)3. 
Later studies of electron-beam and gamma- ray irradi- 



ated solids by Oka and collaborators (Oka 1993 Chan, 
Okumura and Oka|2000 Momose et al|2niiT ) emphasised 
the use ot para-i±2 (i-e. J = 0) enriched solids to achieve 
intrinsically very narrow spectral lines; again the cations 
were interpreted as H;^ and its solvated forms. 

2.2. Broad spectral features 

In parallel with the investigations of radiation-induced 
absorption lines there were studies of broad-band absorp- 
tion in irradiated solid hydrogens. Again t hese studies 



commenced with solids con taining tritium (Richardson 
et al 1981 Poll et al 1983), and later utilised external 



sources of ionising radiation (Selen et al 1984 Forrest 
et al||1990 ). Broad infrared /optical absorptions were re 



ported and interpreted in t erms of bound-bound tran- 



sitions of electron bubbles ( Richardson et al 1981 
et al||1983[ [Selen et al||1984| . And a continuum U 



Poll 

ab- 
sorption was reported and int er preted as the boun d- free 
transitions of electron bubbles ( Forrest et al||199 "o|). 

Surprisingly, irradiated solid hydrogens were also ob- 
served to emit optical/near-inf rared continua in a broad 
band peaking around 825 nm dFor r est et"all|19 92', 'Mag-I 
notta, Collins and Mapoles 1994 F^rrestand Brooks 
1997[ ). Both steady emission and short hashes were seen 
to occur, with essentially the same spectrum. Emission 
from solid deuterium was observed to be much stronger 
than that from solid H2. It was also observed that the 
optical emission increased in response to UV exposure, 
indicating a close connection with the UV c ontinuum ab- 
sorption process ( Magnotta, Collins and Mapoles„1994| ). 
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But the actual emission process was not securely iden- 
tified. One possibility is that the emission is due to ra- 
di ative a.ttachment of e~ and hydr ogen atoms ( Forrest et 
all|l992 [Forrest and Brooks||l997[ ). In gas-phase the H~ 
Bound- tree opacity (e.g. Doughty, Fraser and McEachran 
1966) has a spectral peak in a similar location to the ob- 
served emission from solid D2, but the H~ opacity peak 
is much broader than that of the observed emission. 

2.3. Paramagnetic signals 

By virtue of having no unpaired electrons II2 molecules 
themselves have no Electron Spin Resonance (ESR) sig- 
nal. But the irradiated solid exhibits several features 
(e.g. Kumada et al 2008). First there is a strong signal 
from atomic hydrogen, produced by dissociation of the 
molecules. Secondly there is a signal from free electrons 
in the solid, produced by ionisation of the molecules. 
And thirdly there is a qu artet of lines (Miyazaki, Ya- 
mamoto and Aratono|1995 ) , initially suggested to be due 
to ( Symons||1995), but later recognised as the signa- 
ture of CSymons and Woolley"2000'; 'Suter, Engels 
and Lunell 2001; Kumada, 'I'achikawa and 'I'akayanagi 



2005 p. Further ESR investigations of irradiated solid hy- 



drogens conclusively demonstrated that Hg is indeed the 
source of the quartet lines by studying solid mixtures of 
H2 and HD or D2, for which a variety of isotopomers 
of H;i" arise, each with its own unique ESR signature 
(Kumada, Takayanagi and Kumagai 2006 Kumagai et 

l[lm7\ . — 



2.4. Molecular hydrogen ion clusters 

All of the literature cited in §§2.1,2.2 implicitly as- 
sumes that molecular hydrogen ion chemistry in the con- 
densed phase resembles that of the gas phase, for which 
it is well established that clusters H+ with odd-n oc- 
cur in_jnu£h_greate£_abundance than those with even- 
n ([Clampitt and Gowla nd,1969 Hiraoka 1987 Ekinci, 
Knuth and 'lbennies|2006[ ). I'he siabsequent discovery ot 
Hq in the ESR experiments on irradiated solid hydro- 
gens (§2.3) challenges that assumption. There seems to 
be no experimental data which point specifically to the 
presence of , or its solvated forms, in irradiated solid 
hydrogens: it is possible that the dominant cations are 
Hq isotopomers and their solvated forms. 

The existence of even-n clusters was predicted on the 
basis of ab initio calculations which showed them to 
have comparable binding energies to the odd-n clus- 
ters (Wright and Borkman 1982). The even-n clus- 
ters were then discovered in gas pha se via mass spec- 
troscopy (Kirchner and Bowers 1987), albeit at abun- 



dances very much lower than those of odd-n. Hj was 
seen in these experiments to have much higher abun- 
dance than , or H^g and this result prompted fur- 
ther theoretical investigations of and the other even- 



n clusters ([M ontgom ery and Michels 1987[ [K urosaki 
and Takayanagi 1998ajb|. In particular KmosaET^nH' 



Takayanagi^ (1998ai ) confirmed the stability of and 
identified its most stable configuration as symmetry 
— a structure which is not based on the H;^ moiety, as 
can be seen in figure 1. 

Naively can be thought of as a central with 
an II2 molecule bonded perpendicularly at each end, the 




Figure 1. The D2d configuration of Hg", which is the structure 
corresponding to the global minimum of potential energy for this 
molecule (Kurosaki and Takayanagi, 1998b^ . 

axes of the three moieties being mutually orthogonal. In 
fact each atom in the molecule carries a similar fraction 



of the total charge (Kurosaki and Takayanagi 1998b). 
These authors also demonstrated the stability of larger 
clusters, up to IIj'4. The structure of the clusters was 
shown to correspond to solvation of by H2, with 
the solvating molecules being loosely bound and carrying 
very little charge, and the central moiety only slightly 
perturbed in its charge distribution and bond lengths. 

The existence of very large eve n-n clusters h as been 
demonstrated experimentally by Jaksch et al| p008): 
molecular hydrogen ion clusters as large as n = 120 were 
observed via mass spectroscopy of ionic species formed 
inside ultra-cold (0.4 K) helium nanodroplets. Odd-n 
species were also produced and th e ratio of eyen:odd 
clusters was observed to be small. Jaksch et al (20081 
argued that their high yield of even-n clusters indicated 
that cluster ionisation dynamics are affected by a helium 
matrix. The ionisation dynamics of hydrogen clus ters 
have been studied theoretically by Tachikawa (2000). 



2.5. Summary of impurities in solid H2 

The ESR studies described in §2.3 demonstrate explic- 
itly that irradiated solid hydrogens contain free electrons, 
atomic hydrogen and H^. Free electrons and H-atoms 
can combine to yield H~, and this reaction is strongly 
favoured by the large ratio of the electron affinity of 
hydrogen (755 meV) to the temperature (T < 14 K, so 
kT < 1 meV) . In a low temperature bath of II2 , solva- 
tion of this anion and the cation are both expected. 
Icosahe dral solvation shells have been suggested for bot h 



species (Wang and Andrews 
But it is not clear whether t 



2004 Jaksch et al 2008 



lis remains a preferred gc 



ometry for solvation in the context of a solid II2 matrix, 
where the interface between the solvation shells and the 
H2 lattice presumably plays a significant role in deter- 
mining the form of solvation. 

In the astrophysical context hydrogen is not expected 
to occur as a pure element, but mixed with helium and 
trace amounts of metals. It appears that helium doe s not 
alloy with solid hydrogen (Leventhal and Mills 1991 ), but 
it is possible that a solid hydrogen matrix might contain 
metal impurities. In this paper we ignore that possi- 
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bility and confine our attention to the hydrogen alone 
as there is plenty to occupy us even in this restricted 
case. We further assume that the deuterium abundance 
is Xd ~ 10~^ by number, i.e. comparable to the deu- 
terium abundances measured in diffuse Galactic gas (e.g. 
Sembach 2010). We expect this to manifest itself as an 
abundance of HD : H2 that is Xd by number. 

A solid hydrogen of this composition, when subject to 
ionising radiation, will usually yield that is free of 
deuterium at the time of formation. However, the Zero- 
Point Vibration Energy (ZPVE) of singly-deuterated H^i" 
(plus H2) is about 13 meV lower than that of the un- 
deuterated isotopomer (plus HD) - see Appendix A - 
so at low temperatures condensation of HD onto the H^ 
proceeds. And the net ZPVE is further lowered by a 
similar amount for each of two subsequent deuterations, 
with the result that at temperatures less than about 12 K 
condensation proceeds all the way to (HD);^. 

We use the shorthand (HD);^ to denote the spe- 
cific isotopomer [HD(HD)HD] + , where the central- 
and two end-moieties are encoded in the notation as 
[End(Centre)End] + . 

If the ionisation rate in the solid is sufficiently high 
that the concentration of H^ is of order Xd/3 ~ 3x 10~^, 
then there is insufhcient HD in the lattice to fully deuter- 
ate the H^j" and at these concentrations the singly- and 
doubly-deuterated isotopomers listed in Appendix A will 
constitute a significant fraction of the H^. At still 
higher ionisation rates even these lightly-deuterated iso- 
topomers are in the minority and the H^ isotopomer 
dominates. For simplicity we consider here only very 
high and very low ionisation rates, for which the Hg 



is essentially all in the form of the H, 
topomers, respectively. 



Q and 



HD 



The first two HD condensations on H^ were observed 



by Kumagai et al| (j2007|) in ESR studies of irradiated 
solid hydrogens at 4 K, But these authors did not report a 
detection of (HD);^. As each of the deuteration reactions 
involves lowering the ZPVE by a similar amount there is 
an expectation that if condensation starts it will go all 
the way to (HD) j". We therefore assume that (HD)3 was 
not detected by [Kumagai et al (2007) because its ESR 
spectrum is very complex and/or the lines are broad. 
We note that Kumagai et al (2007) did not report any 
isotopomer of HJ containing D2 or D^ moieties except 
when the solid itself was made from a mixture including 
D2. This demonstrates that the condensation process 
is one where diatomic hydrogen molecules are preserved 
as units which are substituted in or out of the H^ — 
a situation which is unsurprising considering the large 
energy (4.5 eV) that is required to break up H2. 

3. THEORETICAL MODEL OF H+ 

Before tackling H^ itself we performed preliminary 
theoretical studies on H2 and H^. We used various ab 
initio electronic structure met hods, all availabl e within 
the Q-Chcm software package (Shao et al 2006|), to cal- 
culate the bond length and fundamental harmonic fre- 
quency for each molecule. From these calculations it 
was determined that the most suitable level of theory for 
our calculations is CCSD (Coupled Cluster with Single 
and Double excitations), combined with cc-pVTZ basis 
functions (correlation-consistent polarised Valence-only 



Triple Zeta- function). We found that the CCSD/cc- 
pVTZ combination gave a good compromise between ac- 
curacy and computational cost, for both H2 and H^, 
with deviations of less than O.OOlA and 5 cm~^ from our 
most accurate method, CCSD/cc-pVQZ, which utilises 
a Quadruple Zeta-function basis set. Henceforth all the 
results we quote in this paper were obtained by using the 
CCSD/cc-pVTZ level of theory, unless stated otherwise. 

To determine the importance of augmenting our basis 
function set with additional diffuse functions, we per- 
formed calculations of the energy required to break H^i" 
into 2H2+H^. The binding energy at the CCSD/cc- 
pVTZ level is 235.57 kJmol'^ whereas the CCSD/aug- 
cc-pVTZ energy is 236.35 kJmol"^. The effect of aug- 
menting our chosen basis is thus one third of one percent, 
from which we conclude that diffuse functions are not re- 
quired. 

After our calculations on H^ were completed we be- 
came aware of a comprehensive ab initio study of the 
properties of H^ ( Hao et al||2009 ) , which had been pub- 
lished shortly after we conipl eted our initial s urvey of 
the literature. The study of Hao et al (2009) took in 



six different possible structures for H^and utilised var- 
ious levels of theory: both Restricted and Unrestricted 
Hartree-Fock theory, and CCSD with extension to per- 
turbative and full triple-excitations. They used the basis 
sets cc-pVTZ and cc-pVQZ and their augmented coun- 
terparts aug-cc-pVTZ and aug-cc-pVQZ. In all 26 dif- 
ferent combinations of basis-function sets and levels of 



theory were utilised by Hao et al (2009). 

In jj3.2 we compare our results with those of Hao et 



[al] (2009). However, in this paper our emphasis differs 
greatly from theirs in that we have concerned ourselves 
with treating the large anharmonic corrections to the vi- 
bration frequencies of H^, and our study takes in iso- 
topic variations for the reasons given earlier. We under- 
took detailed modelling of only the lowest energy config- 
uration of the molecule, namely the D2d structure. We 
have, however, modelled the D2h transition structure at 
the harmonic level in order to gauge the effects of rota- 
tional delocalisation of the end-groups of the molecule 
(see §3.2.1). 

3.1. Electronic excitations 

Excited electronic states of H^ were investigated us- 
ing the Equations Of Motion formalism (Krylov 2008). 
In addition to the CCSD/cc-pVTZ level of theory we 
utilised CCSD/aug-cc-pVTZ (i.e. augmented by addi- 
tional diffuse basis functions) for this study. At the cc- 
pVTZ level we found the lowest two excited states to 
be 22^1 at 7.44 eV, and I'^E at 8.79 eV. With the aug- 
mented basis these levels shifted to 7.43 eV and 8.35 eV, 
respectively. The state is reached by a valence tran- 
sition and is unaffected by the additional diffuse func- 
tions in the aug-cc-pVTZ basis. The degenerate l^E 
state is stabilized by 0.44 eV in the larger aug-cc-pVTZ 
basis, indicating that it is more diffuse. Not surprisingly, 
given that the binding energy of the ground state of H^ 
is only 2.1 eV (Kurosaki and Takayanagi 1998b), the po- 
tential energy surfaces ot both these states are repulsive, 
so the molecule would dissociate if excited into either. 

We have also investigated the influence of solvation on 
electronic excitations by computing the energies of the 
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two lowest-lying states for Hj[ (112)4. We found them 
to be substantially lower at 6.52 and 7.01 eV (EOM- 
CCSD/cc-pVTZ), but still repulsive. In practice we ex- 
pect much h igher levels of solvation. For example, | Jaksch| 
et al (2008) found evidence for closed icosahedral sol- 
vatfon shells - the first closed shell corresponding to 
H(|'(H2)i2, and the second to (112)12(112)42 - when 
studying very cold (0.4 K) solvated clusters. We there- 
fore expect that the lowest-lying electronic excitations of 
the actual solvated form will be much lower than those 
of 11(1" (112)4 and we assume that some of these states will 
he bound states. This assumption, and the change in the 
electronic spectrum, on solvation, is important in the 
context of excitation of the vibrational modes of (see 
§3.2.2). 

Because we are dealing with a cation, we expect the en- 
ergy required for any further ionisation to be large. We 
have confirmed this by direct calculation: the energy re- 
quired for photo-ionisation ( "vertical ionisation" ) of bare 
H+ is 19.2 eV, and that of H+(H2)4 is 18.3 eV. 

3.2. Vibrational model 

When computing the vibrational properties of a 
molecule the harmonic approximation is usually em- 
ployed. Such calculations are routine with the Q- 
Chem software package and yield both the frequencies 
of the normal modes and the electric-dipole intensities 
of I Aw I — 1 vibrational transitions. In the harmonic ap- 
proximation the Potential Energy Surface (PES) is taken 
to be quadratic about the eq uilibrium geometrv. But 
Hr is a weakly bound system ( Kurosaki and Takayanagi 
1998b), for which the higher-order terms m the PES are 
expected to be relatively large. And furthermore the 
molecule is composed entirely of hydrogens, which are the 
lightest nuclei and thus have the largest vibration ampli- 
tude in any given potential, so the higher-order terms in 
the PES are of particular importance. 

Our model of the effective PES, which we denote UcS, 
is a Taylor expansion in the normal mode coordinates, 
Qi, of the system, around the potential minimum. For 
= 6 atoms there are ZN — 6 = 12 such coordinates 
(the other six degrees of freedom are taken up by rigid 
body rotations and translations) . Our model includes the 
usual quadratic terms, plus contributions that are cubic 
and quartic in the normal mode coordinates. In terms 
of the potential derivatives, r/ (the number of subscripts 
indicating the order of the derivative), we thus have 

^ 3^-6 ^ 3^-6 



Table 1 

Normal modes of the D2d structure of and (HD)j^ 



ijk 



^ 3N-6 

^ 51 VijkiqrqjQkqi- 

ijkl 



(1) 



At the CCSD level only analytic gradients are available, 
so all the 77 terms in equation M had to be computed 
using finite difference methods to evaluate curvatures and 
higher-order derivatives of UcS- We used a step size, h, 
of 0.1 bohr and accurate finite difference expressions that 
give truncation errors of 0(/i^). 

Anharmonic corrections to the vibration frequencies 
were computed using the Vibrational Configuration In- 
teraction (VCI) method (Whitehead and Handy 1975 



Species 


Mode 


Symmetry 


lo 




V 




# 




(km mol^-'^) 


(cm^^) 


(cm^-"-) 


H+ 


12 


Al 


0.0 


3858.1 


3596.0 




11 


B2 


953.3 


3779.1 


3347.6 




10 


Al 


0.0 


2087.4 


1816.0 




9* 


E 


6.9 


1189.5 


1104.8 




8* 


E 


6.9 


1189.4 


1104.7 




7 


B2 


2418.6 


1018.1 


1042.1 




6 


Al 


0.0 


906.9 


884.9 




5* 


E 


13.6 


726.0 


795.0 




4* 


E 


13.6 


726.0 


794.6 




3* 


E 


0.3 


356.2 


794.6 




2* 


E 


0.3 


356.2 


795.0 




1* 


Bi 


0.0 


99.2 


201.3 


(HD)+ 


12 




6.4 


3343.3 


3070.2 




11 




714.3 


3275.3 


2947.7 




10 




274.7 


1827.3 


1574.2 




9* 




5.4 


1073.1 


962.9 




8* 




47.6 


999.1 


908.6 




7 




1142.1 


865.1 


832.5 




6 




299.2 


800.7 


738.1 




5* 




72.9 


616.8 


626.6 




4* 




25.0 


557.4 


613.8 




3* 




2.9 


297.8 


225.3 




2* 




1.2 


291.0 


184.6 




1* 




4.4 


84.7 


120.5 



Note. — Wavcnumbcrs arc given botli with [y') and without {u^o) 
anharmonic corrections. Changes in rotational state of the molecule shift 
the lines away from these wavcnumbcrs. Modes marked with an asterisk 
are strongly affected by rotation of the end-groups, so these modes are 
not good approximations to the eigenmodes of the real molecule (§3.2.1). 



Lin, Gilbert and Gill 2008) in which up to five quanta 
were included. l*'or computational reasons only two-mode 
coupling was considered. That is r\iij terms were consid- 
ered, but 77ijfc terms, where i ^ j ^ k, were not. We also 
excluded all couplings to vibration mode #1 because, as 
discussed later (§3.2.1), this mode is better characterised 
as an internal rotation than a vibration. 

Table 1 shows our results for the normal vibrational 
modes of Hg" . There we list the harmonic frequency, the 
electric-dipole intensity (derived in the harmonic approx- 
imation), and the anharmonic frequency for each mode. 
The four intensities listed as zero are, in the harmonic 
approximation, precisely zero in consequence of the sym- 
metry of the mode (Ai or Bi). Of the remaining eight 
modes, all six of E symmetry are fairly quiet while the 
two B2 modes are both strong. These two modes acquire 
their intensity as a result of displacing the central 
moiety along its axis. Indeed mode 7 can be thought of 
as the central moiety bouncing back and forth along 
the main axis of the molecule and this mode is by far the 
most intense. Modes 12 and 11 both involve stretching 
of the two H2 end-moieties and axial displacement of the 
central atoms. The two modes differ in that mode 12 is a 
symmetric manifestation of these motions whereas mode 
11 is antisymmetric. Consequently the latter involves a 
displacement of the central moiety and thus acquires 
a large intensity, whereas mode 12 simply stretches this 
moiety. Figure 2 shows the motions involved in each of 
the normal modes of 'Hq . 

Comparing the frequencies in the two rightmost 
columns of table 1 shows clearly the high degree of anhar- 
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Figure 2. The normal modes of vibration of Hg . The left column 
shows the even-numbered modes, in descending order starting from 
mode #12 at the top. The right column shows the odd-numbered 
modes, in descending order starting from mode #11 at the top. 
Spheres show the locations of the atoms at rest and red cones 
show the amplitude and direction of motion for each atom. 

monicity in this molecule, exceeding 400 cm^^ for some 
modes. 

Having obtained the PES model as per equation ([T| 
most of the computational work associated with the vi- 
brational calculations is done and it is straightforward to 
investigate the isotopic variants of H^. In §2.5 we noted 
that condensation of HD onto is expected to proceed; 
so in table 1 we also present the normal modes of (HD)j^. 
For this isotopomer the vibrational symmetries of are 
broken and thus there are no silent modes. Indeed modes 
10 and 6 acquire large intensities. Because the effective 
oscillating mass is greater for (HD);^ than for HJ, each of 
the normal modes of (HD);^ has a lower frequency than 
the corresponding mode for H^i" . And whereas Hji" shows 
only two strong lines (in the harmonic approximation), 
there are at least four strong lines for (IID)j]'. The lines 
for both isotopomers are shown in figure 3. 

The normal modes of vibration are not the whole story. 



In the harmonic approximation the normal modes are 
the only modes and they are non-interacting. The inten- 
sities are the same for all |Au| = 1 transitions of a given 
mode, and zero for |Aw| > 2. But we have seen that the 
anharmonicity is large for and we therefore expect 
that there may be significant intensities associated with 
some of the overtones (|Aw| > 2) and combination bands 
(which result from coupling between the normal modes). 
Although we have calculated the frequency shifts due to 
anharmonicity, we are unable to compute the intensities 
of overtones or combinations. 

3.2.1. Errors and uncertainties 

There are various sources of error in our vibrational 
model of H^; although we cannot remove these errors 
it is useful to know roughly how large they are. The 
easiest check to make is whether or not the predicted 
harmonic properties of vary significantly with the 
level of electronic structure theory which is utilise d. That 
tes t can be made by examining the results of |Hao et| 
[al (20091, either in isolation or by comparison wit h our 
own result s. The highest levels of theory utilised by Hao 
([2009) are cc-pVTZ UCCSDT and aug-cc-pVQZ 
and at these levels the harmonic frequencies 



et al 

UCUSb(T; 

typically differ by amounts (mean absolute difference) 
of order llcm~^ and 17cm~^, respectively, compared 
to our own results. These differences constitute a small 
fractional error for most of the vibrational modes. 

Our level of theory corresponds preci sely to the cc- 
pVTZ UCCSD case of |Hao et al] ( |2009| . In this case 
the mean absolute frequency difference between the two 
sets of results is only 6cm~-'^ (or 4cm~-'^ if mode #1 is 
excluded from the compa rison). We attribute these dif- 
ferences to the fact that Hao et al ( j2009) use analytic 



derivatives whereas we employ numerical derivatives. 

In respect of the predicted (non-zero) intensities of the 
various modes there is, however, a substantial di fference 
between our results and those of Hao et al ( 2009 1 . Com- 



pa ring our results t o the highest levels of theory employed 
by Hao et al ( 2009 ) we see good agreement only for mode 
:^11 (diltcrence less than 10%), whereas for mode #7 
their predicted intensity is less than half of ours and for 
modes #8 and 9 it is larger by a factor of 10. The dis- 
agreement is larger still for the four remaining modes, 
their intensities being more than 20 times larger for #4 
and 5, and more than one thousand times larger for #2 
and 3. However, these discrepancies parallel those which 
are evident amongst the various levels of theory of |Hao| 



et al (2009) for each of the modes, to wit: discrepancies 
ot order 10% for mode #11; a range spanning a factor 
of 3 for modes #7, 8 and 9; a factor of 30 for modes #4 
and 5; and a factor of 100 for modes #2 and 3. 
Giv en the huge variation in intensities reported by[ 



et al (2009), for the same mode with different levels of 
theory, we have no cause to doubt our own estimates of 
the mode intensities. 

By comparing our calculations to observed vibrational 
characteristics we are also able to gauge errors beyond 
just the harmonic level. These checks cannot be made 
against itself, because the vibrational frequencies of 
that molecule have not yet been measured. But we have 
checked the accuracy of the quartic approximation in 
equation ([!]) for the case of the H2 molecule: we find 
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Figure 3. Stick diagram of tlie mid-IR lines of Hg" (upper panel), 

and (HD)^ (lower panel). The line strengths are normalised to 
that of mode #7 for each molecule. The E-type modes of the 
structure are not good approximations to the actual eigenmodes 
of the molecules (see §3.2.1) and they are therefore plotted with 
dotted lines. 

that the computed vibration frequency is 61 cm^^ above 
the measured value for that molecule. Employing a sex- 
tic approximation to the PES reduces this error to only 
2cm~^, so the error associated with the quartic approxi- 
mation of equation ([l]) appears to be roughly 60 cm~^ (at 
4,159 cm~^). Another major source of error in our cal- 
culations of vibrational frequencies stems from our use 
of only two-mode coupling in the VCI method. By mod- 
elling Hg we find that the errors introduced by restricting 
our calculations to two-mode coupling could be as large 
as 50cm"i(at 2,521 cm-^). 




Figure 4. The D2h transition structure of Hi" differs from the D2d 
structure of figure 1 only in respect of the dihedral angle between 
the end-groups. This structure is a saddle p oint of the PES and it 
lies only 12 cm~ ^ above the D-^d minimum ( [Kakizaki, Takayanagi| 
[and Shiga|2007l l. 

All of our results here pertain to the isolated 
molecule, and its isotopomers, whereas we expect that 
inside a matrix of solid H2 would be in the solvated 
form, which will have slightly different vibrational prop- 
erties. In §3.4 we give quantitative estimates of the ef- 
fects of solvation on the vibration frequencies of by 
noting the shifts which occur at the harmonic level be- 
tween Hg and H^(H2)4. 

Finally we draw attention to a fundamental source of 
uncertainty in our modelling of the vibrations of H(t" . Al- 
though the structure is the minimum energy config- 
uration of H(l", the structure, in which the two end- 
groups lie in the same plane (see figure 4), is a transition 
structure whi ch lies only 1.5 meV (12cm~^) abo ve the 
minimum (Kakizaki, Takayanagi and Shiga 2007). This 
transition structure can be reached simply by counter- 
rotation of the end-groups of H^ , i.e. the degree of free- 
dom manifest in vibration mode ^\ (see figure 2). But 
the vibration frequency for mode #1 is larger than the 
barrier height and even the zero-point vibration energy 
is sufficient to completely delocalise the structure in this 
coordinate. The resulting probability distribution func- 



tions have b een computed by [Kakizaki, Takayanagi and 



Shiga ( 2007 1 for the isotopomers H^ and and in both 
cases they are almost uniform in the dihedral angle. The 
situation is more complicated in the case of (HD);^, be- 
cause the centre-of-mass of the HD end-groups does not 
coincide with the mid-point of the bond, around which 
the groups rotate. For small amplitude vibrations (in 
mode #1) the result is simply that the H atom oscillates 
with twice the amplitude of the D atom. But that can- 
not be the case for large amplitude motions: if one atom 
performs a complete rotation, so must the other. Rather 
we anticipate that the atoms of the central HD+ group 
will move in response to the rotation of the e nd-groups. 



The (HD)j" isotopome r was not modelled by [Kakizaki 



Takayanagi and Shiga (2007) but we assume that, hke 



Hj and D^, the end groups of the molecule are delo 
calised in the dihedral angle. 

Delocalisation affects the vibrational characteristics of 
the molecule. Strictly speaking it means that none of our 
analysis is valid, because J7efi is a sinusoidal function of 
gi which, over its full range, is very poorly represented by 
equation [ij And to the extent that the properties of a vi- 
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Table 2 

Vibrational frequency shifts for D2h structure 
and solvated form of D2d- 



Mode 


Symmetry 




A, 


# 




(cm~^) 


(cm^-^) 


12 


Ai 


6 


— 50 


11 


B2 


3 


-68 


10 


Al 


10 


-23 


9 


E 


55 


5 


8 


E 


-115 


5 


7 


B2 


-10 


25 


6 


Al 


-3 


27 


5 


E 


120 


34 


4 


E 


-112 


34 


3 


E 


39 


56 


2 


E 


-17 


56 


1 


Bi 







Note. — The symmetry designation refers to 
the isotopomer in the D2d structure. Frequen- 
cies appropriate to the D2h structure are given by 
adding A^h to the values in Table 1. Frequencies 
appropriate to the solvated species arc given by 
adding to the values in Tabic 1. 

bration mode depend on the dihedral angle that mode is 
necessarily coupled to mode #1. All couplings involving 
mode #1 are explicitly excluded from our model because 
mode #1 has more the character of an internal rotation 
than a vibration. 

To gauge the size of these effects we can compute the 
harmonic vibration frequencies for the D2h structure and 
compare them to those of the structure. We have 
done this and in table 2 we show our results as the quan- 
tity Ad,, = UJo{D2h) - UJo{D2d)- 

For mode #1 the D2h structure has an imaginary vi- 
bration frequency because it is a maximum of the PES 
traversed along the coordinate qi. For modes #6,7,10,11 
and 12 we find that the frequencies are the same as those 
of the D2d structure to within 10 cm~^ and the mean ab- 
solute difference is only 6.4 cm^^. These shifts are very 
small and that tells us that those modes are good approx- 
imations to the eigenmodes of the actual system, with its 
rotationally delocalised end-groups. We can therefore be 
confident that we have a reliable estimate of the frequen- 
cies of the corresponding lines. 

But for the remaining six modes, all of which have 
E-type symmetry, the shifts are an order of magnitude 
greater. Degeneracy of the three pairs of E-modcs of 
in the D2d structure is lifted in D2h with splittings 
of 56 cm -1 (modes #2,3), 232 cm^^ (modes #4,5) and 
170cm^^(modes #8,9). Thus the E-modes are strongly 
influenced by the nearly-free rotation of the end-groups 
of the molecule. 

The reason for this striking difference between the E- 
modes and the modes of A and B symmetry is that the E- 
modes all involve displacements of the central two atoms 
perpendicular to the z-axis (i.e. the long axis of the 
molecule), whereas the other modes only displace these 
atoms along the z-axis (see figure 2). Thus E-mode vi- 
brations are affected by the orientation of the end-groups, 
but the other modes are not. 

The large frequency shifts of the E-type modes between 
D2d and D2h structures indicate that the corresponding 
vibrational eigenstates which we have characterised for 



the D2d structure are poor approximations of the ac- 
tual eigenstates of the rotationally delocalised molecule. 
We therefore place correspondingly little weight on these 
modes when comparing our results to observational data. 

3.2.2. Excitation of vibrations 

In solid H2 the temperature is necessarily below the 
triple-point of hydrogen so we have kT < 10cm~^, 
whereas the smallest vibrational quantum of exceeds 
100 cm^^, thus all of its excited vibrational states lie too 
far above the ground state to have any significant ther- 
mal population. There will thus be no thermal emission 
of vibrational quanta from in solid H2 . But the vi- 
brational modes of can be excited by non-thermal 
interactions. The most important non-thermal process is 
likely to be absorption of ambient starlight, for which the 
energ y density peaks at a wavelength of approximately 



^ gt 

1 ^m ( |Mathis, Mezger and Panagia|1983 ), corresponding 



to photon energies ~ 1 eV. 

There are two mechanisms by which may absorb 
ambient starlight: via electronic excitation and through 
vibrational combination bands. Ordinarily one would 
expect the former to be by far the more important of 
these two processes, as combination bands are usually 
very weak. However, there are three aspects of which 
suggest that combination bands may not be completely 
negligible. First is that the energy of the lowest excited 
electronic state is very high at 7.44 eV, where the energy 
density of the ambient radiation fie ld is only one fifth to 
one tenth of i ts value at the peak (Mathis, Mezger and 
Panagia 1983). By contrast there are many vibrational 



combination bands of distributed right across the 
broad spectral peak of the ambient starlight. Secondly, 
even the lowest excited electronic state lies well above the 
vibrational continuum of the electronic ground state of 
, which has a bindi ng energy of only 49 kcal mol~ ^ , o r 



2.1 eV per molecule (Kurosaki and Takayanagi 1998b) 
This means that the bulk of the photon energy cannot 
go into exciting the vibrational modes of . Indeed we 
expect that there will be essentially no excitation of vi- 
brations as the two lowest excited electronic states are 
repulsive and dissociation should occur promptly. Thus 
electronic excitation is an inefficient way of powering vi- 
brational emission from . This is in stark contrast 
with excitation through vibrational combination bands 
— a process which is 100% efficient. Finally we note 
that is an unusually fioppy, anharmonic molecule and 
hence we expect it to have unusually strong combination 
bands. 

Solvated forms of are also fioppy and thus are also 
expected to have unusually strong combination bands. 
However, we have found (§3.1) that solvation to H^(H2)4 
lowers the minimum electronic excitation energy signif- 
icantly. And for greater degrees of solvation we expect 
the energy to decrease further; we assume that this gives 
rise to excited electronic states which are bound. Lower- 
ing the energy of the excited state increases the efficiency 
of excitation of the vibrational degrees of freedom and it 
brings the electronic excitations into the spectral peak 
of the ambient starlight. We therefore anticipate that 
for highly solvated forms of H^, which are prevalent at 
low temperatures, electronic excitation is the dominant 
process. 
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Unfortunately we cannot sensibly model the electronic 
excitation of vibrations, because the excited electronic 
states are sensitive to the form of solvation of the 
molecule. Consequently we are currently unable to model 
the optical absorption spectrum of or the vibrational 
level population subsequent to such absorption. It is the 
latter which determines the strength of any vibrational 
emission because the cascade of emitted quanta will be 
dominated by Av = — 1 transitions in the various modes 
which have been excited. Thus we are unable to model 
the intensities of the mid-IR emission lines from . 

We now consider the band structure which arises from 
rotational transitions of the molecule during emission of 
a vibrational quantum. 

3.3. Rotational hand structure 

In the rigid-body approximation the rotational proper- 
ties of 11(1" are characterised by its moments of inertia or, 
equivalently, its rotational constants. For these are 
B^,By = 1.38 cm-i and = 26.8 cm'^ and for (HD)+ 



= 0.92 cm"^ and B, = 18.5 cm" 



Thus both 



isotopomers resemble a prolate symmetric top. Because 
Bx,By <^ kT < Bz in both cases, rotation about the x 
and y axes may be thermally excited, whereas there will 
be little thermal excitation of rotation around z. 

As emphasised in §3.2.1, is actually a floppy 

molecule and rigid-body rotation results are therefore not 
a good guide to its rotational properties. 

Solvated species of will also be floppy because the 
itself is only slightly perturbed by solvation and the 
bonds between the solvent and solvated molecules are 
themselves weak. The impediment to calculation there- 
fore exists also in this case. A further uncertainty for the 
solvated species is whether or not such a large complex 
is able to rotate at all inside a matrix of solid hydrogen. 
Alternatively, internal rotation of the H^l" might occur 
relative to the solvating shells. 

In view of the difficulties just described we are not able 
to undertake sensible calculations of the rotational band 
structure for . 

3.4. Model spectra 

As we are unable to model either the strength of the 
vibrational emission or the band structure appropriate 
to the various lines, the only real information contained 
in our model spectra is the line frequencies. We therefore 
present the results of our modelling simply by marking 
the wavelengths of the vibration modes in figure 5. Re- 
sults are shown for both and (HD)^. 

We remind readers that modes of E symmetry are 
not good approximations to the eigenmodes of the real 
molecule; they are therefore shown with dotted lines, to 
draw attention to the unreliable nature of the predic- 
tion. We further emphasise the effect of the rotation of 
the end-groups by plotting each mode as a box whose 
spectral width spans the frequencies predicted for the 
D2d and D2h structures (see table 2). For the five well 
characterised modes these boxes are sufhciently narrow 
that they resemble vertical lines in figure 5. 

3.4.1. Effects of solvation 



solvated molecule Hj'4 = H^i" (112)4 in D2d symmetry. Be- 
cause of the much larger number of atoms, 11^4 inevitably 
has many more modes of vibration — 24 more than 
in fact. Of these the highest four frequencies are as- 
sociated with stretching of the four solvent molecules, 
and most of the lowest frequencies are associated with 
vibrations which involve the weak bonds between these 
molecules and the core H^. In between there are modes 
which can be readily identified as perturbed modes #2- 
12 of the core H^. Comparing the frequencies of these 
modes with the harmonic frequencies for bare (table 
1) gives us an indication of the likely shifts associated 
with solvation. Our results are given as the quantity Ag 
in table 2 — a frequency shift which should be added to 
the predictions for the bare molecule in order to obtain 
the vibration frequencies of the solvated form. 

We now apply these shifts to our anharmonic fre- 
quency predictions for and (IID)J^ in order to esti- 
mate the vibration frequencies for the solvated forms of 
these molecules. The results of this procedure are shown 
in figure 5. 



In addition to our theoretical predictions, we have 
show n in figure 5 the IS O mid-IR spectrum of the Orion 
Bar ( Peeters et al||2004 ). This spectrum was chosen be- 
cause tHe'dilarareorEigh quality and the band strengths 
and profiles are representative of the sort of spectrum 
which is most commonly seen, i.e. the "Type A" spec- 
trum jGeballe||1 997( lTokunaga| p97{ [Peeters et ai||2002 
van Diedenhoven ct al 2004). 



I^'or clarity ot presentation we have eliminated from this 
plot sixteen narrow lines arising from atomic or molec- 
ular hy drogen, neon, sulphur or argon. The lines are 
([van Disho eck 2004): H, Brackett-a; four lines from the 
Ptund series ot H; five pure rotational lines from II2; 
[Ar II] 7.0 Aim; [Ar III] 9.0 fiui; [S IV] 10.5 ^im; [Ne II] 
12.8 Atm; [Ne III] 15.6 /xm; and [S III] 18.7 ^im. We thus 
concentrate on the broad bands of observed emission 



(Gillett, Forrest and Merrill 1973]) ; these are often re- 
terred to as the Unidentihed IntraRed (UIR) bands and 
we use that term henceforth. 
It has been proposed that the UIR bands are vi bra- 



tions of Polycyclic Aro matic Hydrocarbons (PAHs) ( Du- 
ley and WiUiams 1981 Leger and Puget|[l984 Allaman- 



Following Kurosaki and Takayanagi 
characterised the harmonic vibration : 



( 1998b I we have 



'requencies of the 



dola et al|1985| ) and this suggestion has been examined in 
considerable detail — see the review by [Tielens (2008). 
However, the identification with PAHs is not completely 
secure. In part that is because there is no predictive the- 
ory for the shapes or relative strengths of the observed 
bands. This means that if the emission from H^ lies 
within these bands then it will be difhcult to place tight 
constraints on the abundance of the isotopomer in ques- 
tion. 

First we compare the Orion spectrum with our results 
for the H^ isotopomer, restricting attention to the five 
lines whose frequencies we can predict with confidence. 
Overall we see no clear correlation between the data and 
our predictions for H^, in either bare or solvated forms. 
For bare H^ mode #6 lies right on the strong 11.3 /.tm 
(885 cm"-*^) band, and for the solvated molecule this mode 
is less than 30 cm~^ from the band, so within the possible 
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Figure 5. I nfrared Space Observatory spectrum of the Orion Bar 
( [Peeters et al 2004). A background fiux oc A^-^ has been sub- 
tracted trom the data so that the broad emission bands can be 
seen more easily. Also plotted here, in the lower half of the figure, 
are the locations of the predicted fundamental vibration lines for 
and (HD)j|' isotopomers in bare and solvated forms. Each line 
is stretched into a band whose spectral width spans the predictions 
for D2d and structures in the case of the bare isotopomers. 

Reliable predictions are shown with solid lines whereas the E-type 
modes are indicated with dotted lines. Predictions are coded thus: 
A ^ H+; B ^ H+(H2)4; C ^ (HD)+; and D ^ (HD)+(H2)4. 

errors (see §3.2.1) both species yield a coincidence. It is 
more difficult to claim a coincidence with any of the other 
modes, either for solvated or bare H^. Unfortunately 
the data shown here do not extend shortward of 3 /im; 
however, we are not aware of any strong emission bands 
in this region. The closest feature to modes #11 and 12 
is thus the 3.3 /.tm emission band seen in the plot and 
this is so far removed from the predicted lines that it is 
plainly not attributable to the isotopomer in either 
bare or solvated forms. Similarly, mode #10 is almost 
200 cm" ^ away from the nearest strong emission band at 
6.2 /im. 

Translating the lack of spectral coincidences into limits 
on the Hq abundance requires information which we do 
not yet have. However, it is fair to say that the data in 
figure 5 give no reason to think that the isotopomer is 
abundant in interstellar space. Although the data shown 
in that figure are an example of a "Type A" spectrum, 
we would have reached the same conclusion if we had 
compared our results to a Type B or Type C spectrum as 
these differ from Type A mainly in respect of the relative 
strengths and sh apes of the various bands, rather than in 
th eir pla cement ( Peeters et al||2002 van Diedenhoven et 
al 2004). It remains possible that spectral coincidences 
could be found in other sorts of data — see the discussion 
of absorption spectra in §5.3. 

The situation is quite different for (HD)|^: figure 5 
shows a striking coincidence between the data and our 
predictions, especially for the solvated form. Coinci- 
dences are seen in the case of the emission bands at 
3.3/im (mode #12), 3.4 /xm (mode #11), 6.2/xm (mode 



#10), 11.3 /im (mode #7), and 12.7 /im (mode #6). And 
we emphasise that this accounts for all of the reliable 
mode predictions for (HD)i^. There is no predicted mode 
coincident with the strong emission bands at 7.7 and 
8.6 /im, but this is of no immediate significance since 
the remaining vibration modes, the E-type modes, are 
not good approximations to the eigenmodes of the actual 
(IID)j]' molecule (§3.2.1). Notwithstanding the uncertain 
identification of the 7.7 and 8.6 /im emission bands, the 
excellent agreement between the data and the well char- 
acterised modes of (HD);^ prompts us to suggest that this 
molecule is responsible for the ubiquitous UIR bands. 



5. DISCUSSION 

Having suggested that (HD);^ is the carrier of the UIR 
bands we now need to consider the broader implications 
of that suggestion. For example, the UIR bands typi- 
cally carry a significant fraction (~ 10%) of the total in- 
tensity seen_on various lines of sight through our Galaxy 



(Onaka 2004), so their carrier must be contributing sub- 
stantiaily to the extinction on these lines of sight. In 
this section we consider (HD)J^ in the broader picture of 
ISM photo-processes. Readers who are unfamiliar with 
current models of interstellar dust m ay wish to consult 



Draine 



Uraine 



(| 2003|) and the many papers in Witt, Clayton and 
. The following discussion is couched en- 
tirely in terms of (HD);^, even though the same physics 
would in most cases apply to the isotopomer. 

5.1. Excitation mechanism 

In §3.2.2 we argued that excitation of the vibrational 
modes of (HD)j|' is expected to occur as a result of ab- 
sorption of ambient starlight. The two contributing pro- 
cesses were identified as electronic excitation and absorp- 
tion in combination bands, with electronic excitation ex- 
pected to dominate for highly solvated forms. Some in- 
sight into the excitation process is available from the data 
in figure 5 where we see that there are peaks at both 3.3 
and 3.4 /im with the former being by far the stronger 
of the two. Our model of (HD)^ does indeed predict 
lines at approximately these locations - modes #12 and 
11, respectively - but the ratio of calculated transition 
strengths is the opposite of that observed: mode #11 
is more than one hundred times stronger than mode 
#12. This suggests that if the observed bands arise from 
(HD)jj' then the excitation process is probably electronic. 
The reason is that the strength of a combination band 
in part reflects the strength of the individual transitions 
involved in that combination, so if these modes were ex- 
cited by absorption in combination bands then we would 
expect that mode #11 would be the stronger of the two 
lines. By contrast, if the excitation is electronic then 
the population of the excited vibrational levels is deter- 
mined by the Franck-Condon factors and is unrelated to 
the transition strengths of the modes. It seems probable 
that the Franck-Condon factors would favour excitation 
of mode #12. The reason is that mode #12 is a symmet- 
ric vibration mode, whereas mode #11 is antisymmetric. 
The Franck-Condon factors reflect the distortion in the 
molecular PES which occurs as a result of the electronic 
excitation, and for the low-lying electronic excitations 
those distortions will be predominantly symmetric. 
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5.2. Relationship to FIR 

In this paper we have concentrated on the vibration 
modes of (IID)j^ itself. But the solvated form has in ad- 
dition a large number of vibration modes which we have 
thus far ignored. For example, in the case of icosahedral 
solvation with a single shell, i.e. (HD)j]'(H2)i2, we would 
expect a total of 84 modes, whereas (HD)j^ has only 12. 
Extrapolati ng from the behaviour seen for l ow degrees 
of solvation ( Kurosaki and Takayanagi|1998b I we expect 
that a complex with a single icosahedral shell would ex- 
hibit 12 high-frequency modes, close to the vibration fre- 
quency of unperturbed H2, and 60 low-frequency modes, 
in addition to the modes common to bare (HD);^. If sol- 
vation extends to a second complete icosahedral shell, i.e. 
(HD);^ (112)12(112)42, then we expect there to be 270 low- 
frequency vibration modes attributable to the solvation. 

Following electronic excitation of the solvated complex, 
courtesy of an absorbed photon, it is probable that many 
of the low-frequency modes will be excited because the 
lowest-lying electronic excitations arc likely to be spread 
over the whole complex. Thus if (HD)j^ is responsible 
for the UIR bands then we expect that Far InfraRed 
(FIR) emission will accompany them. Because of the 
large number of modes crammed into a relatively small 
emission bandwidth, the FIR emission may appear as a 
quasi-continuum if the degree of solvation is high. 

FIR emission is indeed observed to accompany UIR 
band emission and in fact there is generally observed to 
be far more power in the FIR than in the UIR bands (e.g. 
Onaka 2004). In a model which seeks to attribute the 
observed FIR power to solvated (HD);^ that would imply 
that the majority of the observed interstellar extinction 
must be attributed to absorption by (IID)j]". 

5.3. Absorption spectra 

If (HD);^ is present in interstellar space then it should 
manifest itself in absorption as well as emission. The 
relative strengths of the infrared vibrational lines may 
differ greatly between these two processes. Absorption- 
line strengths just reflect the oscillator strengths of the 
transitions (assuming the molecule is in the vibrational 
ground state). These mode intensities are given in table 
1 and plotted in figure 3. Thus we see, for example, 
that in absorption the 3.4 /im line of (HD);^ is expected 
to be much stronger than the 3.3 fim line, whereas in 
emission the 3.3 fim band is observed to be more intense. 
In practice it is unusual to see any of the UIR bands 
in absorption, so there isn't much data to compare our 
results with. However, there is sometimes an aborption 
band seen at 3.4 ju m with no detecta ble 3.3 fim absorption 



strong absorption lines expected for that molecule are 
at 9.6 /xm and 3.0 /im — coincident with known bands 
of strong interstellar absorption. These bands are usu- 
ally attribut ed to amorphous silicates and to water ice, 
respectively fDrainc 2003 van Dishoeck 2004) . It is be- 



accompanying it (Pendleton 2004| . Conventionally this 
band is interpreted in terms ot aliphatics — long-chain 
hydrocarbon molecules, very different to the polycyclic 
species which are usually invoked to explain the other 
UIR bands. 

We are not aware of any other interstellar absorption 
bands coincident with our mode predictions for (IID)j^. 
However, as our predictions match the locations of the 
observed UIR emission bands it seems unlikely that this 
situation presents any particular difficulty for interpret- 
ing the UIR emission in terms of (HD)|^. 

Looking again at figure 3, but this time at the pre- 
dictions for the isotopomer, we see that the only 



yond the scope of this paper to construct an astrophysical 
model in which these absorption bands are interpreted 
in terms of H^ ; here we simply note the spectral coinci- 
dences. 

We have already noted that the optical absorption 
spectrum of solvated H^ is probably dominated by elec- 
tronic excitations (§3.2.2) and that this expectation is 
reinforced by the data in figure 5 (§5.1). We are unable 
to calculate the electronic absorption spectrum, but we 
can give a qualitative description as follows. 

The large number of vibrational modes possessed by a 
single solvated molecule means that we can expect to see 
a very large number of absorption lines. Taking the bare 
molecule, for example, there are 12 modes of vibration. 
Assuming the Franck-Condon factors favour only two dif- 
ferent quantum states for each mode gives us 2^^ ^ 4,000 
different vibrational states and the same number of ab- 
sorption lines. For the solvated species this number is, 
of course, vastly greater: for a single icosahedral shell, 
i.e. (HD);[(H2)i2, it is ~ lO^^. Given that the binding 
energy of the clu sters is only about 2.2 eV in the gr ound 



electronic state (Kurosaki and Takayanagi 



1998b), the 



spacing between individual lines is expected to be typi- 
cally very small so that they create, in effect, a continuum 
absorption. 

Individual absorption lines should nevertheless be de- 
tectable because the large number of vibration modes 
also leads to a large dynamic range in line strength. 
For example, if we assume, as before, that the Franck- 
Condon factors favour only two different quantum states 
in each mode, but that one of these states has typically 
twice the Franck-Condon factor of the other, then the 
intensity ratio between the strongest and weakest lines is 
of order the total number of lines. 

All of the foregoing arguments apply equally well to 
any bound, excited electronic states. It is likely that 
there will be a number of such states and this corre- 
spondingly multiplies the number of possible absorption 
lines. 

An unusual aspect of (HD);^ is its high levels of an- 
harmonicity. One result of this is that in any vibrational 
sequence the level spacings are not uniform. Thus, even 
in cases where the Franck-Condon factor is large over a 
broad range of one vibrational quantum number, there 
will be no regular spacing in the corresponding sequence 
of absorption lines. 

If we try to make a connection to the astronomical data 
then we are looking for correlated optical continuum and 
line absorption phenomena in which the optical lines are 
unidentified. We therefore consider the possibility of a 
connection to the Diffuse Interstellar Ban ds (DIBs). The 
properties of DIBs have been reviewed by [ Herbig ( 1995 1 
and Sarre (2006), for example. Their phenomenology is 
rich and we do not attempt a detailed assessment vis- 
a-vis (HD)j^ absorption; instead we confine ourselves to 
the following points. First, DIBs are known to corre- 
late well with optical continuum absorption, as gauged 
by reddening. Secondly, DIBs are concentrated in the 
optical band, consistent with a carrier that has a small 
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binding energy. Thirdly, although there have been sug- 
gestions of regularities in DIB spacings (notably, Herbig 
1988), with the best available data it appears that there 
are no statistically significant regulari ties in line spacing 
for intervals smaller than 400 cm" ~i (Hobbs ct al 2008) 
20091. This is consistent with absorption by a hoppy 
molecule. 

Finally, the lack of very strong correlations amongst 
the vario us DIB strengths, measured on different lines 
of sight ( [McCall et "al||2010[ [Friedman et al||2011| , is a 
puzzling teature ot the data. But it can be understood 
in the context of an absorption spectrum in which any 
individual absorption line represents a miniscule fraction 
of all the absorptions, as we expect for solvated (HD)j]". 
This is because that fraction will depend on the smallest 
details of the wavefunctions of the ground- and excited- 
states and the dependencies will be different for each 
absorption line. Branching ratios may therefore exhibit 
a strong environmental dependence. 

We also note that DIB strength is known to correlate 
with the colum n-density of atomic hydrogen, not m olec- 
ular hydrogen (Herbig 1993 Friedman et al 2011). At 
first sight this might seem to argue against (HD)j]" be- 
ing the carrier, given its intimate association with solid 
molecular hydrogen. However, the correlations are deter- 
mined using measurements that are sensitive to diffuse 
gas and insensitive to precipitates. Furthermore it is pos- 
sible that the diffuse atomic hydrogen could actually arise 
from photodissociation of H 2 which is itself undetected 
( [Allen and Diaz-Miller||2004[ ). 

Finally we recall that (HD)|^ resides in a matrix of 
solid H2, which contributes absorption, as do the H^ 
ions which we expect are also present (§2.5). The elec- 
tronic transitions which make up the Lyman and Werner 
bands of H2 are very strong and they are expected to in- 
troduce high levels of absorption in the Far UltraViolet 
(FUV) . This is qualitatively consistent with the observed 
rapid rise in extinction toward the FUV (e.g. Fitzpatrick 
2004) . Although the properties of the Lyman and Werner 
transitions haye been w ell studied in the gas phase of H2 
( Abgrall et al|[1993a|b ), we are not aware of any studies 
ot their appearance m the condensed phase. Similarly, 
we are not aware of any studies of the properties of H~ 
in a solid hydrogen matrix. Quantitative estimates of 
their contributions must therefore be deferred. 

5.4. Red fluorescence and the ERE 

In §2.2 we noted that irradiated solid hydrogens are 
seen to fluoresce in the far red and near infrared. This flu- 
orescence may be of interest in an astrophysical context 
as a possible explanation for the Extended Red Emission 
(ERE) which has been observed f rom various types of 



nebulae and the broader ISM — see [Witt and Vijh| ( [2004t 
for a review of the ERE phenomenon. With re ference to 
the r e quirem ents on the ERE carrier listed by [Witt and[ 
Vijh (2004) we note the following points of similarity 
between the two phenomena: (i) the spectral peaks oc- 
cur in similar places; (ii) the spectra are similarly broad, 
unstructured and confined to similar wavelength ranges; 
(iii) both are photoluminescence phenomena; and (iv) 
the ERE arises from an abundant element with Galaxy- 
wide distribution. There is no doubt that hydrogen is 
an abundant element with a Galaxy-wide distribution. 



We do not yet know if it is found in solid form in astro- 
physical contexts, but if it is then it would be a strong 
candidate for the carrier of the ERE. 

6. CONCLUSIONS 

It is remarkable that a single, small hydrogen molecule 
exhibits most of the strong mid-IR bands which per- 
vade astronomical environments. This fact alone makes 
(HD)j]" of interest to astronomers, even though we can- 
not yet place it in our theoretical picture of the interstel- 
lar medium with any degree of confidence. The 7.7 and 
8.6 /im astronomical bands, which are often the strongest 
of all the observed UIR features, are absent from our 
predictions for (HD);^. But there are several vibrational 
modes which we are unable to characterise, leaving open 
the prospect that some of those modes are responsible 
for the bands in question. 

Inside a matrix of solid hydrogen, (HD)^ is solvated 
and the many weak bonds in the resulting cluster are 
expected to give rise to a large number of FIR vibra- 
tional transitions. Like the mid-IR bands, these cluster 
modes will be excited following electronic absorption of 
optical photons, and we expect a correspondingly close 
relationship between mid-IR and FIR emissions from the 
solvated (HD)3" molecule. 

We arc unable to compute the electronic absorption 
spectrum of solvated (HD);^ because it is sensitive to 
the details of the solvation, but the floppy nature of 
the molecule leads us to expect line properties which 
are broadly consistent with those observed for the DIBs. 
And the multitude of weak lines yields an apparent con- 
tinuum of optical absorption. 

To these phenomena one can add FUV absorption con- 
tributed by the solid hydrogen matrix itself and a red 
fluorescence, possibly attributable to H~. Thus, in an 
astronomical context, solid hydrogen is expected to man- 
ifest itself in a variety of ways, each of which seems to 
resemble one of the observed astronomical phenomena 
collectively attributed to "dust." We conclude that seri- 
ous consideration should be given to the possibility that 
solid H2 is abundant in the interstellar medium. 
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Table 3 

Harmonic and anharmonic Zero-Point Vibration Energies 
for various Hg" isotopomers 



Isotopomer Harmonic ZPVE Anharmonic ZPVE 

(mcV) (mcV) 



[H2(H2)H2] + 


1010.0 


1000.9 


[H2(H2)HD] + 


964.5 


955.3 


[H2(HD)H2] + 


962.3 


952.7 


[H2(HD)HD] + 


916.6 


906.2 


[H2(DH)HD] + 


915.9 


907.5 


[HD(H2)HD] + 


918.6 


909.5 


[HD(HD)HD] + 


869.9 


860.7 


[H2(H2)D2] + 


913.6 


905.8 


[H2(D2)H2] + 


912.3 


903.2 


[H2(D2)D2] + 


813.9 


807.3 


[D2(H2)D2] + 


815.9 


810.2 


[D2{D2)D2] + 


714.4 


709.8 



Note. — A nhar monic corrections to the ZPVE were computed 
using equation | A 1 1 with Z^i^^ and (j^ijk terms neglected. 

APPENDIX 
ISOTOPIC VARIANTS OF H+ 

In §2 we noted that the dominant isotopomer of H^i" in astrophysical environments is usually either or (HD)j^. 
In the laboratory, however, conditions may be chosen which favour different isotopomers. Here we tabulate properties 
of the ten isotopomers which, in addition to and (HD)!^, arise in solids made from mixtures of II2 with HD or D2. 

In table Al we list the Zero-Point Vibration Energies of all t hese isotopomers. Thes e values for the ZPVE include 
anharmonic corrections, treated as a perturbation according to ( Schuurman et al||2005 ) 



37V-6 



ZPVE = - 
2 



1 

32 



3iV-6 

E 

ijk 



<t'iik4'kjj 



1 

48 



3^-6 
i]k 



^ijk 



1 

32 



3N-6 



z, 



kin 1 



(Al) 



(A2) 



where the various ^'s are the mode-coupling constants, given by 

Vijk 

and similarly for the quartic constants. We have not evaluated the three-mode coupling coefficients {4>ijk, see §3.2) and 
so we ignore the third term on the right-hand-side of equation (Al). We also ignore the last term (Zkm), associated 
with rotation, as here we are treating mode #1 as a vibration on a par with the other modes. 

To determine whether or not isotopi c condensation o ccurs it is also necessary to know the ZPVE for H2 and HD or 
D2. From table 39 in the appendix of |Herzberg| ( |1950[ ) we find ZPVE(H2)=268.8 meV, ZPVE(HD)=233.7 meV and 
ZPVE(D2)=191.4 meV. 

In tables A2 and A3 we list the wavenumbers and intensities of the normal vibrational modes of ten isotopomers, 
to complement the results for and (HD);^ given in table 1. Likewise the plots in figures Al and A2 complement 
those in figure 3. The theoretical modelling is as described in §3. We remind readers that the line intensities reported 
here are computed in the harmonic approximation, even though we have calculated anharmonic corrections for the 
frequencies of the lines. All transitions in tables A2 and A3 are w = 1 O 0, for the mode in question, and rotation is 
neglected. The notation for an isotopomer consisting of Endi, Centre and End2 moieties is [Endi(Centre)End2]^. 
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Table 4 

Normal modes of the 0-24 structure of Hg" isotopomers 
containing H2 and HD groups. 



Species 


Mode 




Wo 


V 




# 


(km moI~^) 


(cm~^) 


(cm-i) 


[H2(H2)HD] + 


12 


422.7 


3821.0 


3530.2 




11 


395.6 


33U7.y 


2973. D 




10 


0.6 


2081.4 


1774.8 




9* 


8.2 


1187.6 


1076.4 




8* 


13.7 


1169.0 


1074.3 




7 


2282.6 


1004.6 


1000.3 




6 


69.2 


860.7 


828.7 




5* 


11.4 


723.7 


740.8 




4* 


40.3 


627.3 


738.4 




3* 


1.5 


349.0 


227.5 




2* 


6.2 


332.6 


740.8 




1* 


2.2 


92.9 


170.6 


[H2(HD)H2] + 


12 


0.0 


3858.0 


3566.2 




11 


992.2 


3777.8 


3346.5 




10 


274.6 


1842.6 


1573.2 




9* 


3.6 


1104.6 


1024.1 




8* 


19.7 


1020.9 


970.0 




7 


797.1 


927.3 


937.5 




6 


810.5 


862.1 


830.0 




5* 


14.9 


725.7 


724.6 




4* 


24.8 


660.7 


724.6 




3* 


2.9 


325.0 


727.5 




2* 


3.0 


320.2 


211.2 




1* 


0.0 


99.2 


179.2 


[H2(HD)HD]+ 


12 


431.3 


3820.9 


3541.4 




11 


426.1 


3306.2 


2970.2 




10 


249.8 


1832.0 


1624.8 




9* 


2.8 


1074.7 


1001.6 




8* 


22.5 


1018.3 


925.4 




7 


1056.9 


923.7 


855.7 




6 


509.2 


805.4 


789.5 




5* 


12.4 


723.5 


664.3 




4* 


56.1 


566.8 


666.8 




3* 


1.0 


313.4 


202.9 




2* 


3.2 


307.5 


175.4 




1* 


1.8 


92.8 


145.8 


[H2(DH)HD]+ 


12 


445.4 


3819.9 


3558.5 




11 


Ann '7 
409.7 


3307.5 


2979.5 




10 


301.1 


1837.6 


1611.3 




9* 


4.3 


1103.3 


990.7 




8* 


37.6 


1002.3 


964.4 




7 


1446.6 


866.9 


926.9 




6 


29.7 


857.9 


778.5 




5* 


21.6 


658.1 


694.7 




4* 


64.4 


615.9 


694.7 




3* 


0.8 


314.6 


675.7 




2* 


3.5 


298.3 


270.4 




1* 


1.4 


92.7 


163.0 


fHDrHTlHDl+ 


12 


6.4 


3343.6 


3123.5 




11 


677.2 


3276.9 


2942.5 




10 


0.1 


2075.5 


1748.5 




9* 


3.5 


1167.6 


1037.9 




8* 


29.3 


1166.3 


1041.3 




7 


2293.7 


979.4 


941.7 




6 


2.5 


829.3 


775.8 




5* 


57.8 


643.0 


684.1 




4* 


2.5 


600.6 


677.0 




3* 


10.7 


328.3 


217.3 




2* 


4.0 


322.4 


187.6 




1* 


5.6 


84.9 


133.5 



Note. — Wavenumbers are given both with (i/) and without 

(cuo) anharmonic corrections. Changes in rotational state of the 
molecule shift the lines away from these wavenumbers. Modes 
marked with an asterisk are strongly affected by rotation of the 
end-groups, so these modes arc not good approximations to the 
eigenmodes of the real molecule (§3.2.1). 
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Table 5 

Normal modes of the D2c( structure of Hg" isotopomers 
containing H2 and D2 groups. 



Species 


Mode 


lo 




V 




# 


(km mol^^) 


(cm~^) 


(cm"'^) 


[H2(H2)D2] + 


12 


465.5 


3819.2 


3476.3 




11 


205.2 


'-7r\ 

27U3.8 


2474.5 




10 


2.6 


2075.0 


1785.0 




9* 


9.5 


1185.8 


1064.7 




8* 


0.9 


1153.5 


1052.5 




7 


2126.6 


996.9 


979.4 




6 


250.3 


760.0 


749.0 




5* 


9.5 


721.4 


708.7 




4* 


2.2 


581.4 


708.7 




3* 


4.2 


341.2 


699.7 




2* 


2.2 


310.5 


248.2 




1* 


0.0 


88.2 


149.4 


[H2(D2)H2] + 


12 


0.0 


3857.8 


3483.6 




11 


1039.1 


3776.2 


3350.8 




10 


0.0 


1521.4 


1361.6 




9* 


17.6 


923.5 


874.8 




8* 


17.6 


923.4 


874.8 




7 


1331.9 


832.0 


773.9 




6 


0.0 


880.3 


880.7 




5* 


16.1 


648.1 


665.0 




4* 


16.1 


648.0 


664.6 




3* 


9.4 


303.4 


182.9 




2* 


9.4 


303.4 


184.0 




1* 


0.0 


99.2 


154.4 


[H2(D2)D2] + 


12 


504.5 


3817.7 


3506.9 


11 


252.9 


2700.8 


2470.3 




10 


6.1 


1497.8 


1345.7 




9* 


21.7 


920.8 


850.3 




8* 


1.5 


846.7 


826.1 




7 


682.0 


852.2 


804.4 




6 


587.7 


690.0 


660.2 




5* 


9.8 


640.8 


610.9 




4* 


12.1 


520.3 


583.7 




3* 


1.9 


283.5 


240.9 




2* 


2.0 


270.7 


164.0 




1* 


0.0 


88.2 


122.9 


[D2(H2)D2] + 


12 


0.0 


2730.1 


2600.0 




11 


390.6 


2677.7 


2443.0 




10 


0.0 


2063.6 


1718.6 




9* 


2.1 


1149.0 


1015.2 




8* 


2.1 


1149.0 


1015.2 




7 


2321.7 


927.9 


876.9 




6 


0.0 


648.7 


629.2 




5* 


0.5 


574.9 


626.2 




4* 


0.5 


574.8 


625.9 




3* 


9.6 


297.5 


223.9 




2* 


9.6 


297.4 


224.5 




1* 


0.0 


70.1 


98.5 


fDo('Do")DDl+ 


12 


0.0 


2729.2 


2540.4 




11 


477.0 


2673.3 


2451.7 




10 


0.0 


1476.6 


1311.5 




9* 


3.5 


841.4 


773.7 




8* 


3.5 


841.4 


773.7 




7 


1210.3 


720.2 


687.2 




6 


0.0 


641.6 


605.3 




5* 


6.8 


513.6 


537.9 




4* 


6.8 


513.6 


537.7 




3* 


0.1 


252.0 


213.6 




2* 


0.1 


252.0 


214.0 




1* 


0.0 


70.1 


85.6 



Note. — Wavenumbers are given both with {y) and without 

(cjo) anharmonic corrections. Changes in rotational state of the 
molecule sliil'l the lines away from these wavenumbers. Modes 
marked with an asterisk are strongly affected by rotation of the 
end-groups, so these modes arc not good approximations to the 
eigenmodes of the real molecule (§3.2.1). 



Interstellar solid hydrogen 



17 



1.2 
^ 1.0 


L [H2(H2)HD]^ 


-I 


"D 

I 0.8 


: 


: 


^ 0.6 






(u 0.4 






^ 0.2 


^ I I 




0.0 







8 9 10 




Figure 6. Stick diagram of tlic mid-IR lines of five isotopomers of H^. The line strengths are normalised to that of mode #7 or #11 for 
each molecule. The E-typc modes of the D2d structure axe not good approximations to the actual eigenmodes of the molecules (see §3.2.1) 
and they are therefore plotted with dotted lines. 
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Figure 7. Stick diagram of the mid-IR lines of five isotopomers of Hg . The line strengths are normalised to that of mode #7 for each 
molecule. The E-type modes of the D2d structure are not good approximations to the actual eigenmodes of the molecules (see §3.2.1) and 
they are therefore plotted with dotted lines. 



